(/X , e ) Conversion; A Symbiosis of Particle 
and Nuclear Physics 

T. S. KOSMAS and J. D. VERGADOS 
Theoretical Physics Division, University of loannina, GR 45110, Greece 

Abstract 

(/U~ , e~ ) conversion is the experimentally most interesting lepton 
flavor violating process. From a theoretical point of view it is an 
interesting interplay of particle and nuclear physics. The effective 
transition operator, depending on the gauge model, is in general de- 
scribed in terms of a combination of four terms of transition operators 
(isoscalar and isovector, Fermi-like as well as axial vector-like). The 
experimentally most interesting ground state to ground state transi- 
tion is adequately described in terms of the usual proton and neutron 
form factors. These were computed in both the shell model and RPA. 
Since it is of interest to know the portion of the strength exhausted 
by the coherent (ground state to ground state) transition, the total 
transition rate to all final states must also be computed. This was 
done i) in RPA by explicitly summing over all final states ii) in the 
context of the closure approximation (using shell model and RPA for 
constructing the initial state) and iii) in the context of nuclear matter 
mapped into nuclei via a local density approximation. 

We found that, apart from small local oscillations, the conversion 
rate keeps increasing from light to heavy nuclear elements. We also 
find that the coherent mode is dominant (it exhausts more than 90% 
of the sum rule). Various gauge models are discussed. In general the 
predicted branching ratio is much smaller compared to the present 
experimental limit. 
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1 Introduction 



All currently known experimental data are consistent with the standard 
model (SM) of weak and electromagnetic interactions. Within the framework 
of the SM, baryon and lepton quantum numbers are separately conserved. In 
fact one can associate an additive lepton flavor quantum number with each 
lepton generation which appears to be conserved. There are thus three such 
conserved quantum numbers Lg, and Lj. each one associated with the 
lepton generations (e~, z/g), i/^), (r~, i/^), with their antiparticles having 
opposite lepton flavor. It is in fact these quantum numbers which distinguish 
between the three neutrino species if they are massless. 

Most theorists, however, view the SM not as the ultimate theory of na- 
ture but as a successful low energy approximation. In possible extensions of 
the SM it is legitimate to ask whether lepton flavor conservation still holds. 
In fact in such gauge models (Grand Unified Theories, Supersymmetric Ex- 
tensions of the SM, Superstring Inspired Models) such quantum numbers are 
associated with global (non local) symmetries and their conservation must 
be broken at some level. 

Motivated in part by this belief the search for lepton flavor violation, 
which began almost half a century ago (Hincks and Pontecorvo, 1948 [|^], 
Lagarigue and Peyrou, 1952 0, Lokanathan and Steinberger, 1955 0, see 
also Frankel, 1975 |^) has been revived in recent years and is expected to 
continue in the near future. In the meantime the number of possible reac- 
tions for testing lepton flavor has been increased. The most prominent such 
reactions are 

/i ^ 67 (1) 

r — > 67 and t ^ fi'j (2) 
fi ee~^e~ (3) 
r — > ee'^e~ , r fj,e^e~ (4) 
r 6/i^/i~, r yU/i^yU~ (5) 
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Kl /i^e^, K+ TT+fie (6) 
(/x^e^) ^ (/i^e^) muonium — antimuonium oscillations (7) 

fi~ + [A, Z) —>■ e~ + [A, Z) [muon — electron conversion) (8) 
Finally, one could have both lepton and lepton flavor violating processes 

like 

{A, Z) (A, Z±2) + e^e^ {^^^v - decay) (9) 

/x^ + {A, Z) ^ e'^ + {A, Z — 2) {muon — positron conversion) (10) 

From an experimental point of view the most interesting reactions are 
(|I]), (H), (IP and (0). None of the above processes has yet been seen. 
The best limits obtained are 

R = r(^+ e+7) ^ ^ ^Q_n ^9Q% (^^^ 
r(/i+ -> e+Ve^^) 

set by LAMPF (Bolton et ai, 1988 0), 

- %;ZuP < 1-° X 1°-" (12) 
set by SINDRUM at PSI (Bellgardt et al, 1988 §), 

i?e-^ = ^ 'J'^ < 4-6 X 10-1^ (90% CL) (13) 

1 (/i — >■ a//) 

set by TRIUMF (Ahmad et al., 1987 0) using a Time Projection Counter 
(TPC) and 

R...= ^^'~lr"'''m'^^ < 4.4x10-^^ (90% CL) (14) 
1 (yU ^ all) 
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set by SINDRUM II (Badertscher et ai, 1991 [g). For a Pb target the 
corresponding (/i~,e~) conversion limit is (Ahmad et ai, 1988 0) 

Re-N < 4.9 X 10"^° (90% CL) (15) 

From a theoretical physics point of view the problem of lepton flavor 
non-conservation is connected with the family mixing in the leptonic sector. 
Almost in all models the above process can proceed at the one loop level 
via the neutrino mixing. However, due to the GIM mechanism in the lep- 
tonic sector, the amplitude vanishes in the limit in which the neutrinos are 
massless. In some special cases the GIM mechanism may not be completely 
operative even if one considers the part of the amplitude which is indepen- 
dent of the neutrino mass (Langacker and London, 1988 0, Valle, 1991 [|1^ 



Conzalez-Garcia and Valle, 1992 |n|). Even then, however, the process is 
suppressed if the neutrinos are degenerate. It should be mentioned that 
processes (|l])-(@) cannot distinguish between Dirac and Majorana neutrinos. 
Process (P) can proceed only if the neutrinos are Majorana particles. 

In more elaborate models one may encounter additional mechanisms for 
lepton flavor violation. In Grand Unifled Theories (GUT's) one may have 
additional Higgs scalars which can serve as intermediate particles at the one 
or two loop level leading to processes ([I|)-(|^). In supersymmetric extentions 
of the standard model one may encounter as intermediate particles the su- 
perpartners of the above. Lepton flavor violation can also occur in composite 
models, e.g. technicolor |T^. In fact, such models have already been ruled 



out by the present experimental bounds (see eqs. ([Tl| ) - ([T5|) ). 

The observation of any of the processes (|I|)-(|TO|) will deflnitely signal new 
physics beyond the standard model. It will severely restrict most models. It 
may take, however, even then much more experimental effort to unravel spe- 
ciflc mechanisms responsible for lepton flavor violation or flx the parameters 
of the models. The question of lepton flavor non-conservation has been the 
subject of several review papers (Scheck, 1978 [|1^, Costa and Zwirner, 1986 



Engfer and Walter, 1986 []T3|, Vergados, 1986 |]TB|, Bilenky and Petcov, 



1987 [0, Melese, 1989 p|, Heusch, 1990 |T9|, Herczeg, 1992 Schaaf, 
1993 Kosmas, Leontaris and Vergados, 1994 [^). In the present re- 
view we will focus our attention on recent theoretical developments of the 
subject. We will pay little attention to the experimental situation, since we 
do not intend to dublicate the experimental review which recently appeared 
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(Schaaf, 1993 1^). Furthermore, the reader can find an interesting account 



of the early experiments by Di Leila [2^ . 

From a nuclear physics point of view the most interesting muon number 
violating process is the (/i~,e~) conversion in eq. (|). Its sister (/i~,e"'") 
conversion in eq. (0) is much more complicated, since it involves two nucle- 
ons. Furthermore, lepton violation is more likely to be seen in neutrinoless 
double beta decay of eq. (P). In this report we will focus our attention on the 
e~) conversion. We will concentrate mostly on the evaluation of nuclear 
matrix elements. As we have mentioned above, experimentally the most im- 
portant transition is that to the ground state. It is, however, important to 
know which fraction of the total strength goes to the ground state. We will, 
therefore, also evaluate the total transition rate to all final states in various 
nuclear models. 

The (/i~,e~) conversion, compared to its competitors /i 67, /i ^ 3e 
etc., has some experimental advantages p2| : 



1. The detection of only one particle is sufficient. No coincidence is 
needed. 

2. For electrons with the highest possible energy, i.e. Eg ^ m^c^ — e;, with 
eb the muon binding energy, the reaction is almost background free. Indeed 
the reactions of induced background are: 

i) Muon decay in orbit. There is a tiny tail in the region of interest which 
is proportional to (-E^"^ — E^)^, i.e. very small {E^""^ denotes the maximum 
energy of the background electrons). But in this region the shape is known 
and it can be subtracted out. 

ii) Radiative muon capture. Indeed this can be a source of background 
since the photon can decay to e+e~ pairs as 



Ij- + {A, Z) {A,Z -1) + iy^+ 7 

I 



+ e- (16) 



If the neutrino and the positron carry away zero kinetic energy, the back- 
ground electron can be confused with interesting electron. We notice, how- 
ever, that the maximum electron energy for process ([T6|) is 



EZ"" = m^c' -et- m,c^ - A = E, - A - m.c^ (17) 
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where A is the difference in the binding energy of the two nuclei involved in 
eq. (p^BI). By a judicious choice of the target nucleus A can be quite large 
(e.g. A = 2.5MeV for ^^C). Thus, one has a background free region if one 
restricts oneself to the coherent mode (the final nucleus of eq. (P) is left in 
its ground state). 



2 Mechanisms for lepton flavor violation 

Obviously, in all models which allow ^ ej to proceed reactions 
of eqs. (D), (P) and (p!0D can also proceed via a virtual photon (see fig. 1). 
The dot in these figures indicates that the vertex is not elementary but very 
complicated. Such processes are called photonic. Muon-electron conversion 
can, of course, occur via mechanisms which do not involve the photon (non- 
photonic mechanisms). These are Z-exchange diagrams, i.e. the photon 
of fig. 1(b) is replaced by a Z-boson. In most models these diagrams are 
less important than the photonic ones. Another possibility is an effective 4- 
fermion contact interaction (box diagrams). We notice that here it is possible 
to have both protons and neutrons participating which may lead to a nuclear 
enhancement (see fig. 2). 

The most popular scenario for (/i^,e^) conversion involves intermediate 
neutrinos, see fig. 3 for the photonic case and fig. 4 for the box- diagram 
case. In these diagrams the neutrinos which propagate are the mass eigen- 
states. The neutrinos involved in weak interactions, z/g and u^, are not sta- 
tionary states but linear combinations of neutrino mass eigenstates given by 
U^j and Uej. Both light Uj and heavy Nj intermediate neutrinos can prop- 
agate. Since, however, the matrix U is unitary, the part of the amplitude 
which is independent of the neutrino mass vanishes (GIM mechanism). Thus, 
the leading contribution is proportional to Am^/M^, for light neutrinos, or 
Am1fM^/m%, for heavy neutrinos [0. This means that the process is sup- 
pressed for neutrinos with masses very different from the W-boson mass. 
This is discouraging, since all reasonable models do not yield neutrinos in 
the region of the W-mass. 

Another possibility is to enlarge the Higgs sector. One can induce {fi~ , e~) 
conversion in models with two Higgs isodoublets. The most dominant con- 



tribution, however, is expected to occur at the two loop level [22|. This 



calculation has only been done for the /x ^ 67 reaction. Such calculations 
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show that it is possible to reasonably adjust the parameters of the model so 
that a branching ratio not far from the experimental limit can be obtained. 
We should note, of course, that for real photons one encounters only the elec- 
tric dipole and magnetic dipole form factors (see next section). For (/i~,e~) 
conversion (virtual photons) one also needs the monopole form factors which 
have not been calculated. One expects, however, to be able to obtain large 
branching ratios when this is accomplished for fast yU — > 67 decay. 

Other exotic scalars, like the singly charged S~^, the doubly charged isos- 
inglet S^^ and isotriplet x"*"^; can also lead to flavor violations Such 



enlargements are not, however, favored by recent theoretical models, espe- 
cially those which are superstring inspired, and they are not going to be 
further discussed. 

Another interesting possibility is the supersymmetric extension of the 
standard model. One now has additional particles which are the superpart- 
ners of the known particles. With those extra particles participating as in- 
termediaries one can have a plethora of new diagrams [^, |2^. The most 
important for (/i~,e~) conversion are shown in fig. 5. 

The mixing matrix V entering the vertices is 

V = StSe (18) 

where Se is the unitary mixing matrix for the charged leptons and Sg is the 
corresponding one for the S-leptons (see ref. |^). To leading order (tree 
level) the matrices Se and Sg are the same, i.e. V is diagonal, and lepton 
flavor violation occurs. If, however, one goes beyond the tree level and takes 
into account renormalization effects, the matrix V is no longer diagonal so 



that one can have lepton flavor violation |22 . 



3 Expressions for the amplitude of (/i , e ) 
conversion 



The amplitude for the (yU , e ) conversion can be cast in the form [£2 



47m (1) , C (2) A /-.q^ 
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where the first term is the photonic and the second the non-photonic contri- 
bution, q is the momentum transfer and C, takes the form 



— } W-boson mediated models ^^q'^ 

777,3^2/"^!) Supersymmetric models 

777.3/2 is the gravitino mass and rria is the relevant s-quark (supersymmetric 
partner of quark). 

The hadronic currents are 

= iV7;,i±Z3iV, {photonic) (21) 



4'^ = Nlx\ [(3 + /S/yrs) + {fv + /A/3r3)75] iV, {non - photonic) (22) 
(A^ = Nucleon) while the leptonic currents are 



j'l) = (/mi + l,fEi)i<y^''— + (/eo + 75/Mo) 7" - ^1 (23) 

= ^^(Pi)7' ^ (/y + Ia1,)u{p,) (24) 

where /3 = /3i//9o, is the ratio of the isovector to the isoscalar component of 
the hadronic current at the quark level. The form factors /mi, /bo, /mo, 
/y and as well as the parameter (3 depend on the assumed gauge model 
and the mechanism adopted (the isoscalar parameter (3^ is absorbed in the 
definition of the leptonic form factors). For purely left-handed theories the 
number of independent form factors is reduced in half since 

/bo — —/mo, /bi — —/mi, /a = —fv (25) 
In some models involving the 1^-boson one has 

/eg /ei 



(26) 



while in the supersymmetric models discussed above one finds 
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1 Tfl^ 

Ieo = -/mo = -i^m'^9{x)^r- (27) 

1 171?' 

AnafEi = -A-nafui = --fjo^ f{x)^f- (28) 

fv = -fA = -^m'Mx)^ (29) 

/ "^3/2 



and 



0„ = Ul'4, A = i-i4 (30) 

9 9 mj- 9 9 m j 

d a 

The functions g{x), f{x), fb{x) depend on the quantity 7 = m^/ma- Since, 
however, this quantity is much smaller than unity, we get 

f{x) ~ 9{x) ~ h{x) ^ ^ (31) 

4 Effective nuclear transition operator and 
nuclear matrix elements 

The first step in constructing the effective transition operator is to take the 
non-relativistic limit of the hadronic currents in eqs. (pi]), (^). This leads 



to the operators p2|, p6[] 



A A 

^0 = E (3 + fv^r,,) e-'^-^^ , 1^ = -^^^ + /?r3, ) e"^^-"-^ (32) 



with ^ = fyjfj^ and 



= -, = 0, /y = 1, /3 = 3 (photonic case) (33) 
6 



5'v = = 7:, fv = 1, /a = 1-24 (non — photonic case) (34) 



9 



For neutrino mediated processes 



^ _ 1 30 ^ _ 1 25 light neutrinos 

1 1 ' 15/6 heavy neutrinos 

i.e. (3 = 5/6 ~ 0.8 in both cases. For the supersymmetric model discussed 
above 

/3o^5/9, /3i = ^, I.e. /? = 0.6 (36) 
The factor l/Vs in fl was introduced for convenience. Thus one has 

\ME\' = \< f\Qo\ ^,^i>\' +3/1 \< f \ n \ fi >\^ (37) 



The second step is to factor out the muon Is wave function ||23, ^ i.e. 



|< m^,{r)\z >\'^<K >\<f\n\t >\' (38) 



where 



,2 ^_ /rfV|$M(r)lV(r) 



with $/^(r) the muon wave function and p(r) the nuclear density. The above 
approximation was recently found to underestimate the width in heavy nuclei 
by as much as 40% (see below sect. 5). However, it is still a good approx- 
imation for the branching ratio Re-N of gQ- (0); since it affects both the 
numerator and the denominator by the same way. 

With the above operators one can easily proceed with the evaluation of 
the relevant nuclear matrix elements. We distinquish the following two cases. 

4.1 The coherent (^u , e ) conversion matrix elements 



For 0^ —>■ 0^ transitions only the operator Qo in eq. (37) contributes. One 



finds [22 



<^\no\^>=~gv{S + fvP)F{q^) (40) 

where 
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F{q') = F,{q^) + l-^Fiv(g^) (41) 
J + JvP 

Fz{q') = ^J dhp,{v)e~^^-\ F^iq') = ^J p^{v)e-"'^^-^ (42) 

Fz and Fn are the proton, neutron nuclear form factors with Pp(r), Pn(r) 
the corresponding densities normahzed to Z and N, respectively. 
Then, the branching ratio Re-N takes the form 



1 r I "^M 1 ~ |2 i"^^ 1 ~ |2 T 

(43) 

where k and -^ph carry all the the dependence on the nuclear physics, i.e. 
and 

In eq. ( ^41) is a combination of the coupling constants entering the ordi- 
nary muon capture (G^ ~ 6) and fpR is the well known Primakoff function 
[p7| , |28| which adequately describes the ordinary muon capture throughout 
the periodic table. It is approximately given by 



/pR ~ 1.6|- 0.62 (46) 

It is sometimes convenient to factor out the nuclear dependence from the 
dependence on the rest parameters of the theory ||2^, i.e. we write 

Re-N = PI (47) 

where the quantity p is independent on nuclear physics. The quantity 7 takes 
the form 



11 



^ G'ZfpniA, Z) ^^""^ 
In the case of the supersymmetric model discussed above 7 takes the form 



7 = (3 + Y j ^P^- (49) 

At this point it is of interest to compare the branching ratio of (/x~,e~) 
conversion to that of the yU — > 67 decay. One finds 



+ ^ Ivh (50) 



i?e7 GtT \3 4 

4.2 Total {fi~,e~) conversion matrix elements 

As it was mentioned in the introduction, only the coherent rate is of experi- 
mental interest. It is, however, important to know what portion of the sum 
rule is exhausted by the coherent mode. We thus have to evaluate the total 
matrix element 



ML = fvT. (£) l< / I ^0 I ^ >r +3/1E O l< / I ^ 
where 



(51) 



g/ = - C - {Ef - Egs) (52) 

Ef, Egs are the energies of the final and ground state of the nucleus. This 
evaluation clearly can be done in a model in which the final states can be 
explicitly constructed. This is a formidable task, however, and only in simple 
models can easily be done (e.g. RPA PO, 



The other alternative is to use some approximation scheme. The first is 
the so-called "closure approximation" PT|. In this approximation one 



first replaces the momentum qj by a suitable average, i.e. qj ^ k =< qf >. 
Thus p2|, 



no(q/) ^ E ^o(j>-"^"^- = no{k) (53) 
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17(q;) ^ Y: ^(J>-'^-^^ = ^^(k) (54) 
j 

where 

^oO') = 3 + /yT3,-, = (e + /9T3,)a/x/3 (55) 

Then, we write 

(If 



Sa = l</l^^(q/)l^> 



- 2 E<^l"'"WI/>-</l"(^)M> (56) 

or using closure over the final states we get 

Sa = ^ < i \ ^'^(k)-Q(k) \i >= — \ <i \ Q,ib I i > + < i I O26 I i > [ 
ml m2 1- J 

(57) 

where 

= E^^(i)-^(^') (58) 
j 

^6 = E^^(^)-^(J>""''^'*"'^'^ (59) 
The computation of Sy is analogous with 

^^16 = E ^o''(j>o(j) (60) 
j 

= E a;o+(z)^o(j>-'''-('-^-^^-^ (61) 

We thus find 

Mlt = /^5y + 3/1^^ (62) 
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The real question is what is the appropriate value of k? In earlier calcula- 



tions ||25|, ^ an average excitation energy E = 20MeV was used which was 
derived from the ordinary muon capture phenomenology. It was recently 
realized, however, that the correct excitation energy must be appreciably 



smaller ^6] due to the presence in (/i~,e~) conversion of the coherent mode 
with Ef = 0. This channel is absent in muon capture. It seems reasonable, 
therefore, to use an average energy defined by 

E,(^/|</|nU.>l^ 

This definition cannot in practice be used since presupposes knowlendge of 
the final states, which we like to avoid. It can, therefore, be estimated in a 
simple model like RPA (see next section). In the context of RPA we can also 
estimate easily the effect of the gs ^ gs correlations by using the Thouless 
theorem and defining a correlated vacuum in terms of the uncorrelated one 
(see ref. p2| , p6|). This shows that both the coherent as well as the total 
(/i~, e~) conversion matrix elements are a rescaling of the uncorrelated ones 
and that the ground state correlations strongly reduce all matrix elements 
(see sect. 5). 

Another method for calculating the {fi~,e~) conversion widths has re- 
cently been proposed by the Valencia group This method em- 
ploys nuclear matter mapped into nuclei via a local density approximation 
(L.D.A.) utilizing the relativistic Lindhard function. It has been proved that 
this method reproduces very well the ordinary muon capture data. In this 
process only the incoherent channels are open. The method, however, has 



been also applied to the incoherent part of (/i~,e~) conversion ||3^. The 
Lindhard function was computed taking into account particle-hole excita- 
tions of p-n type in a local Fermi sea without invoking the approximation 
of eq. (pSl). The coherent mode is evaluated independently by using a local 
density approximation in eq. (pTf). 



5 Results and discussion 

As we have mentioned in the introduction, in the present report we focus 
on the nuclear dependence of (/i~,e~) conversion rates. In this section we 
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present and discuss the nuclear matrix elements for both the ground state 
to ground state transition (coherent mode) and the sum over all final states 
(total rate) obtained in the context of the three methods discussed above, 
i.e. (i) shell model, (ii) quasi-particle RPA and (iii) nuclear matter mapped 
into nuclei. 



5.1 Coherent (^u , e ) conversion 

For the coherent mode one essentially needs only calculate the proron and 
neutron nuclear form factors. In table 1 the results for the 0"*" — > 0+ transition 
matrix elements obtained in the framework of the shell model and quasi- 
particle RPA are shown. They describe the mechanism for the photonic 
diagrams {(3 = 3) and non-photonic ones (/3 = 5/6) discussed in sect. 2. 
We see that the two methods give about the same results. In both methods 
we have followed the conventional approach and ignored the muon binding 
energy (see table 1 results labeled shell model and QRPA(i)). To check the 
validity of this approximation the calculation was repeated in the context 
of RPA (see table 1 results labeled QRPA(ii)) by explicitly calculating the 
binding energy e?,. For heavy nuclei it was found that the consideration of 
eft gives about 40% larger results. This is expected, since in this case the 
momentum transfer at which the form factors are calculated is smaller. 

We mention that similar results have also been found by using with the 
method of nuclear matter mapped into nuclei These results are shown 



in table 2. With the latter method the accuracy of using the approximation 
of eq. (^) was also tested by doing exact calculation of the muon-nucleus 
overlap which is possible in the context of the nuclear matter mapped into 



nuclei method |^ . We see that the exact calculations for the absolute rates 
(case II of table 2) differ appreciably from those of the approximation of 
eq. (|38|) (case I of table 2). We should mention, however, that we expect 
this approximation to be much better in the branching ratio provided that 
the total (/i~, Ufj) rate is calculated in the same way. One then can use the 
Primakoff's function (see eqs. (H3) and (ff^)) which was derived under this 
approximation. 

The effect of the ground state to ground state correlations was estimated 
in the context of RPA as was discussed in sect. 4.2. The results obtained 
for '^^Ti nucleus are shown in tables 3 and 4. We see that these correlations 
reduce both the coherent and total matrix elements by about 30% which is 

15 



in agreement with other similar results ||3^, 
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5.2 Incoherent (^u , e ) conversion 

The incoherent {fi~ , e~) conversion involves the matrix elements of all the ex- 
cited states of the participating nucleus. These have been calculated with two 
methods: (i) by explicitly constructing the final states in the context of the 
quasi-particle RPA |26| (see tables 3 and 4 results labeled QRPA(explicit)) 
and (ii) by summing over a continuum of excited states in a local Fermi sea 
in the framework of the nuclear matter mapped into nuclei method (see 
results labeled L.D.A. in table 3). 

The total rates have also been obtained by assuming closure approxima- 
tion and employing a mean excitation energy of the studied nucleus as was 
discussed in sect. 4.2. The results obtained this way both in the context of 



shell model 251 and RPA sum rules 261 are shown in tables 3 and 4 for vari- 



ous mean excitation energies E (see sect. 4.2 for details). For RPA the mean 
excitation energy used was found by firstly calculating one by one the excited 



states included in eq. (^). Since such a calculation in the context of the 
shell model is very tedious, the value of E used in shell model sum rules was 
chosen in analogy to that of the muon capture reaction. The big difference 
appeared between the two values is justified because, in the (/i~,e~) con- 
version the dominant gs gs transition (not existing in the muon capture) 
contributes only in the non energy weighted sum rule of eq. ([53|). The above 
large mean excitation energy used in shell model sum rules overestimated the 
incoherent matrix elements (especially in the region of heavy nuclei), since 
they depend strongly on E. 

As we have stressed above, among all the open channels for {fi~,e~) 
conversion the transition to the ground state is of experimental interest. A 
useful quantity is the ratio of the coherent rate divided by the total (/i~, e~) 
rate, i.e. 



1 ME 


2 

coh 


1 ME 


|2 

\tot 



V = (64) 



By dividing the results obtained for the coherent mode (Mg^^g^) and the 
total (yU~,e~) matrix elements (Mj^^J, either they are given directly from a 
sum rule or by adding independent coherent and incoherent results, we found 
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that the coherent channel dominates (see tables 3 and 4) and that throughout 
the periodic table rj > 90% (see ref. and [0). 



6 Summary and conclusions 

In this report we have investigated the {A, Z){fi~ ,e~){A, Z)* reaction em- 
phasizing its dependence on nuclear physics. Using reasonable assumptions 
it was possible to factor out the nuclear dependence. The nuclear matrix 
elements of great interest were those entering the ground state to ground 
state transition. These were described in terms of proton and neutron form 
factors. Shell model calculations give similar results with RPA [^]. The 
predicted proton form factors agree with those extracted from electron scat- 
tering. As expected |^ ^ ground state correlations reduce the rates by as 
much as 30%. The predicted proton form factors agree with those extracted 
from electron scattering. 

The total matrix element was also computed. First, by explicit calcula- 
tions involving all the final states in the context of RPA. Second, by employ- 
ing a sum rule with a suitable mean excitation energy E ^ 2MeV, i.e. much 
smaller than previously expected p5|]. This was done both in shell model 



and RPA. Again the ground state correlations tend to decrease the predicted 
rates. For the calculation of the total rate we also used a new method which 
employs nuclear matter mapped into nuclei via a local density approxima- 
tion. Needless to say that, though all three methods basically agree with 
each other, the accuracy of such calculations is not as good as those for the 
coherent process. It is, however, pretty certain to expect that throughout 
the periodic table the coherent mode dominates, i.e. rj > 90% This, as we 
have mentioned in the introduction, is experimentally very important. 

Returning to the coherent production we find that there is some depen- 
dence on nuclear physics. In fact we find that 

1.67 <K< 1.89 (65) 

(see eq. (^5|) ) which is a small effect. The variation in 'jph is, however, much 
more pronounced. 

1.6 < %h < 26 (66) 
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The (/i , e ) conversion rate does not show a maximum around A ^ 60, as 



it was previously beheved [§81 , but keeps increasing all the way up to the 
heaviest elements ^ . 

The nuclear dependence on {A, Z) is not strong enough to overcome the 
extra power of a in eq. (0). We find 

1.5x10-3 <^<3xl0-2 (67) 

(the largest refers to ^^"^Sn and the lowest to ^"^C). 

Finally, it is not encouraging that the predicted branching ratios for the 
muon number violating processes are much smaller than experiment. In the 
super symmetric model described above we find 

1.2 X 10-^^ < Re~N < 2.4 X 10-^^ (68) 

We must stress, however, that this should not discourage the experimental- 
ists, since we do not as yet have a complete theory to adequately describe 
such processes. 
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Table 1: Coherent , e ) conversion matrix elements M^^^^^ calculated 
in the context of shell model and quasi-particle RPA (see text). 



Nucleus 


Photonic Mechanism (/5 = 3) 


Non-Photonic Mechanism (/3 = 5/6) 




Shell Model 


QRPA (i) 


QRPA (ii) 


Shell Model 


QRPA (i) 


QRPA (ii) 


22^^26 


142.7 


135.2 


139.6 


374.3 


363.2 


375.2 




187.5 


187.8 


198.7 


499.6 


498.2 


527.4 




212.9 


212.7 


227.8 


595.8 


596.2 


639.5 


112^ J 
48 '-^"64 


274.2 


280.0 


346.7 


769.4 


785.8 


983.3 


70 O92 


313.6 


311.0 


484.3 


796.0 


840.3 


1412.1 


208 pu 
82 -r Oi26 


240.2 


287.5 


582.9 


631.4 


767.5 


1674.9 



Figure Captions 
Fig. 1 

Photonic diagrams for the elementary process /i — > ej, fig. 1(a), and the 
induced processes (/x~,e~) conversion, fig. 1(b), /j, —>■ e+e", fig. 1(c) and 
muonium-antimuonium oscillations, fig. 1(d). • is a complicated vertex. 



Fig. 2 

Box-diagrams of (//~, e~) conversion in an effective 4-fermion contact in- 
teraction. Both protons and neutrons participate. 

Fig. 3 
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Table 2: Coherent widths (in arbitrary units) for the photonic and non- 
photonic mechanisms obtained with the exact muon wave function (column 
labeled I) and with the approximation of eq. (38) (column labeled II). 



Nucleus non- photonic mechanism (/3 = 5/6) photonic mechanism (/3 = 3) 



A 


Z 


I 


II 


I 


II 




12 


6 


0.52 10"^ 


0.51 10-^ 


0.21 10-^ 


0.21 10" 


-4 


24 


12 


0.91 10"^ 


0.90 10-3 


0.38 10-3 


0.37 10" 


-3 


27 


13 


0.14 10-2 


0.14 10-2 


0.56 10-3 


0.55 10" 


-3 


32 


16 


0.29 10-2 


0.28 10-2 


0.12 10-2 


0.12 10" 


-2 


40 


20 


0.64 10-2 


0.62 10-2 


0.27 10-2 


0.26 10" 


-2 


44 


20 


0.72 10-2 


0.69 10-2 


0.26 10-2 


0.25 10" 


-2 


48 


22 


0.94 10-2 


0.91 10-2 


0.35 10-2 


0.34 10" 


-2 


63 


29 


0.21 10-^ 


0.19 10-^ 


0.78 10-2 


0.72 10" 


-2 


90 


40 


0.47 10-^ 


0.42 10-^ 


0.17 10-1 


0.16 10- 


-1 


112 


48 


0.62 10-^ 


0.54 10-^ 


0.22 10-1 


0.19 10- 


-1 


208 


82 


0.13 10° 


0.89 10-1 


0.41 10-1 


0.29 10" 


-1 


238 


92 


0.12 10° 


0.73 10-1 


0.35 10-1 


0.22 10" 


-1 
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Table 3: Total and coherent (// , e ) matrix elements for ^^Ti (non-photonic 
mechanism, P — 5/6). The ratio 77 of eq. (64) is also shown. 



Method 




M^ 

qs—tqs 


E 


ML 


v{%) 


L.D.A. 




375.7 




410.5 


91.5 


Shell Model {sum — rule) 


374.3 


20.0 


468.0 


80.0 


QRPA (explicit) 




363.0 




386.4 


93.9 


QRPA{sum — rule) 




363.0 


0.5 


366.2 


99.1 


QRPA{sum — rule) 




363.0 


5.0 


376.5 


96.4 


QRPA{sum — rule) 




363.0 


20.0 


382.8 


94.8 


QRPA + Corr (sum - 


- rule) 


236.2 


0.5 


238.6 


99.0 


QRPA + Corr {sum - 


- rule) 


236.2 


5.0 


245.4 


96.3 


QRPA + Corr {sum - 


- rule) 


236.2 


20.0 


249.5 


94.7 



Tabic 4: Total and coherent {fi ,e ) matrix elements for the photonic mech- 
anism {P = 3). See caption of table 3. 



Method 




M^ 

qs—^qs 


E 


ML 


v{%) 


Shell Model{sum — rule) 


144.6 


20.0 


188.8 


67.2 


QRPA {explicit) 




135.0 




146.2 


92.3 


QRPA{sum — rule) 




135.0 


1.7 


138.3 


97.6 


QRPA{sum — rule) 




135.0 


5.0 


140.6 


96.0 


QRPA{sum — rule) 




135.0 


20.0 


141.7 


95.3 


QRPA + Corr {sum — 


rule) 


87.8 


1.7 


90.4 


97.1 


QRPA + Corr {sum — 


rule) 


87.8 


5.0 


91.8 


95.6 


QRPA + Corr {sum - 


rule) 


87.8 


20.0 


92.6 


94.8 
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Photonic diagrams oi {/j, ,e ) conversion in a model with mixed inter- 
mediate neutrinos. 

Fig. 4 

Box diagrams (for protons and neutrons) of e~) conversion in a model 
with mixed intermediate neutrinos. 

Fig. 5 

SUSY diagrams for (//~,e~) conversion. Intermediate charged s-lepton 
mixing. 



24 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9408011vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9408011vl 



I^L Vj (Nj) 




photon from charged leptons 




(A,Z) (A,Z-1) (A,Z) (A,Z) (A,Z+1) (A,Z) 



